Two-dimensional Rydberg gases and the quantum hard squares model 
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We study a two-dimensional lattice gas of atoms that are photo-excited to high-lying Rydberg 
states in which they interact via the van-der-Waals interaction. We explore the regime of dominant 
nearest neighbor interaction where this system is intimately connected to a quantum version of 
Baxter's hard squares model. We show that the strongly correlated ground state of the Rydberg 
gas can be analytically described by a projected entangled pair state that constitutes the ground 
state of the quantum hard squares model. This correspondence allows us to identify a first order 
phase boundary where the Rydberg gas undergoes a transition from a disordered (liquid) phase to 
an ordered (solid) phase. 
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The study of quantum many-body models with long- 
range interactions is a rapidly growing and flourishing 
field [IHS] fueled by recent successes in creating ultra cold 
polar molecules [H [5] as well as strongly interacting ultra 
cold gases of Rydberg atoms [3 [7]. Rydberg gases of- 
fer widely tunable inter-atomic interactions which allow 
to explore correlated quantum matter in vastly different 
regimes, ranging from supersolids to exotic spin models 
[SHTS], A particularly interesting regime is achieved at 
densities where the strong interaction between Rydberg 
states competes with the excitation laser. In the extreme 
case which is illustrated in Fig. [ija), the simultaneous 
excitation of atoms located within a spatial radius i?b 
(called the "blockade radius") is forbidden [19]. This 
blockade effect renders Rydberg atoms into "hard ob- 
jects" whose mutual exclusion induces strong correlations 
in the ground state. In a one-dimensional setup this was 
analyzed in Ref. [18j , where it was shown that the ground 
state can even be obtained analytically for certain combi- 
nations of the laser parameters. A current challenge is to 
understand and to identify interesting and useful features 
of strongly correlated two-dimensional quantum systems. 
Since the situation in two-dimensions is substantially 
more involved, exact analytical solution of strongly cor- 
related two-dimensional quantum systems are scarce. 

In this work we present a detailed study of a two- 
dimensional spin lattice gas motivated by recent experi- 
mental achievements in creating Rydberg lattice systems 
|20| and the future perspective to implement quantum 
spin systems with polar molecules Our system is 
intimately related to a quantum version of a hard ob- 
jects model - Baxter's hard squares f2T1 - whose ground 
state represents a coherent spin state with built-in near- 
est neighbor exclusion and belongs to the class of pro- 
jected entangled pair states (PEPS). We show that for 
certain parameters this strongly correlated state provides 
an accurate analytical description of the two-dimensional 
Rydberg gas. This highlights the possibility of creating 
PEPS in a physical system with binary interactions which 
is important as such states are considered a valuable re- 




FIG. 1. (Color online) (a) Atoms arranged on a square lattice 
and coherently photo-excited from their ground state \g} to 
a Rydberg state |e) with Rabi frequency i} and detuning A. 
In the regime of ideal nearest-neighbor blockade a Rydberg 
atom entirely blocks the excitation of its four nearest neigh- 
bors (crosses). Exciting a Rydberg atom is, thus, equivalent 
to placing a hard square (dashed lines) centered on the posi- 
tion of the Rydberg atoms on the lattice. The shaded area 
indicates the range of the Rydberg-Rydberg interaction, (b) 
The phases of the quantum hard squares model [eq. ([2|] for 
different values of the parameter ^ defined in the text. For 
C < Ccrit: liquid phase in which both sub-lattices (solid and 
open circles) are populated. For ^ > ^crit : ordered phase with 
broken sub-lattice symmetry. 



source for quantum information processing [531 US] ■ Fur- 
thermore, our study yields insights into the phase struc- 
ture of the Rydberg gas allowing us to identify a first 
order quantum phase transition from a disordered quan- 
tum liquid to an ordered state [Fig. [l|b)]. 

The system under consideration consists of an ensem- 
ble of atoms regularly arranged on a square lattice and 
coherently photo-excited from their ground to a Ryd- 
berg nS'-state by a laser with Rabi frequency O and de- 
tuning A [Fig. [ija)]. We focus on the situation where 
the system is prepared in a Mott insulating state with 
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one atom per site in a deep large-spacing optical lattice 
[24| . The interaction between ground state atoms at dif- 
ferent sites is negligible, while Rydberg atoms strongly 
interact with each other via a van-der-Waals interaction 
Uij = Ce/\ri — rj\^. Here, r^, rj denote the position vec- 
tors of the excited atoms, Cg is the interaction coefficient 
and we define the interaction strength at nearest neigh- 
bor distance a as V = Cq/q^. In the following, we will 
consider the regime -Rb ^ a shown in Fig. [l^, where the 
interaction is so strong that the simultaneous excitation 
of nearest neighbor particles is completely blocked. As 
demonstrated in [TH] , it is then convenient to work in the 
interaction picture with respect to the nearest neighbor 
interaction, since the Hamiltonian in this picture explic- 
itly displays the effect of the excitation blockade. 

Neglecting small corrections of order fl"^ /V the Hamil- 
tonian reads, 

k.m k,rn 

V s-^ 

) ■ (1) 

The Pauli matrix cr^ ^ changes the state of the lattice 
site {k,m) = Ti/a from down to up and vice versa 
with the identification |t) = |e) and \i) = \g). The 
projector n^^m = 0, 1 determines if the site is occu- 
pied by a Rydberg atom, and we use its complement, 
— i^k,rm to define a plaquette operator 'Pk,m = 

Pk-l,mPk,m~lPk+l,mPk,m+l- 

In the following, we will refer to the terms of H that 
solely depend on the density Uk^m as "classical terms". 
They describe the energy penalty to be paid if the laser 
is detuned by A from resonance, and the interaction be- 
tween Rydberg atoms, where we have cut off the long- 
range tail of Uij at next-nearest neighbors, since its 
strength quickly decreases with increasing distance of 
the particles. In contrast, we will call the first term of 
eq.([l]) "quantum term". It describes the coherent (de- 
)excitation of atoms at a rate J7 and is a five-particle 
interaction term: A particle at site (k, m) can only be 
excited, if all of its four nearest neighbors are simultane- 
ously in the ground state. As indicated in Fig. [TJa), the 
action of the quantum term corresponds to placing a hard 
square on the site of the excitation or removing it from 
that site. Therefore, the realization of a two-dimensional 
Rydberg lattice gas in the ideal nearest neighbor blockade 
regime is equivalent to implementing a quantum lattice 
model of interacting hard squares which are created and 
annihilated by the laser. 

The many-body ground state of Hamiltonian ([T]) is 
not known analytically. Its numerical determination by 
a straightforward diagonalization of H is restricted to 
small system sizes, due to the exponential growth of the 
Hilbert space dimension with increasing particle number. 
We therefore propose an alternative method to study the 



ground state of the strongly interacting Rydberg lattice 
gas. The main idea is to consider an auxiliary model 
system whose many-body ground state wave function is 
known analytically and which possesses the same symme- 
try properties and the same quantum term as H. Using 
this wave function as a variational input state it is then 
possible to approximately determine the phase diagram 
of our original problem. 

Our candidate model system is a quantum version of 
Baxter's classical hard squares model [21; on a square lat- 
tice. Its Hamiltonian is of the so called Rokhsar-Kivelson 
type |25j . i.e. it can be written in the form, 

i?HS = ^^^4,™^fc.™> (2) 
k.rn 

with positivc-scmidefinite, self-adjoint operators. 



hk,m. = y _^ ^ [^k,m + C ^"-fc.m + C(l - "fc,m)] ^fe,m , 

and a non-negative parameter ^. The mathematical 
properties of the hk^m entail that the ground state en- 
ergy of -ffns is zero and, consequently, the ground state 
wavefunction |^) is annihilated by all hk^m- It is given by 

10 = ^ exp I -C ^ alJ^k,n I |0) , (4) 

y'^i \ k,m J 

where |0) denotes the empty lattice (all atoms in the 
ground state), a'^^^^ = i'^k,m + *'^fc.m)/2 is the spin rais- 
ing operator at site (fc, to) and the normalization con- 
stant. Physically, the many-body ground state of 77hs 
is a coherent superposition of all many-particle states 
in which Rydberg atoms are arranged in configurations 
compatible with the nearest neighbor (or hard squares) 
exclusion. Each configuration is weighted by a factor 
(~C)"i where n is the number of occupied lattice sites. 
The state |^) is a PEPS and can be seen as a coherent 
spin state with built-in nearest-neighbor exclusion. 

The normalization constant is the partition func- 
tion of the classical hard-squares model with fugacity 
Hence, the parameter ^ controls the occupation of 
the lattice. Using the partition function, classical ob- 
servables like the density of particles and density-density 
correlations can be calculated. Interestingly, the expec- 
tation value of non-classical observables can also be cal- 
culated using Z^, since they are connected to expectation 
values of classical observables taken in the ground-state 
1^). For example, the mean values of cr^ „ is given by 

(eK™io^KJ? = -2r^(nfe,™)c. 

The classical hard squares model displays an order- 
disorder phase transition between a liquid and a solid 
phase with broken sub-lattice symmetry }21j . The ther- 
mal fluctuations that drive this transition are mapped to 
quantum fluctuations by the PEPS solution Q. There- 
fore, the quantum hard squares model will also undergo 
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FIG. 2. (Color online) Main plot: Energy of the 6 x 6 Ry- 
dberg lattice gas on the ^-manifold (c.f. red dashed line in 
Fig. [3|. Red line: exact ground state energy of (T|. Black 
line: energy calculated using the coherent state (Ml with the 
same parameter ^ as taken for the parameterization of the 
^-manifold. Squares: Energy calculated by taking the PEPS 
with an independent parameter rj that minimizes {rj\H{£,)\ri) . 
Upper inset: Same for the fraction of excited atoms. Lower 
inset: Overlap between the exact ground state of Hamiltonian 
^ with the PEPS Q using the parameter ^ (black) and the 
optimized parameter rj (squares). For details see text. 



a liquid-solid transition at a critical value ^crit ~ 1-9. 
Here however, it is driven by quantum fluctuations [26 
[c.f. Fig. [ijb)]. Since a order-disorder quantum phase 
transition is also anticipated for ultra cold Rydberg gases 
[SHE] and since the Hamiltonians a nd ([2| have the 
same symmetry properties, the state (|4|) appears to be a 
natural candidate for performing a variational study of 
the phase diagram of the Rydberg lattice gas. 

To solidify this statement we will now quantify the in- 
timate connection of the Rydberg lattice gas problem 
with the quantum hard squares model. By multiply- 
ing out the quantum hard squares Hamiltonian ^ and 
comparing it termwise to eq. ([T]), one finds that the 
quantum term is identical in both cases. The classical 
terms of iJns that are linear in nk.m are identical to 
those of H if we choose the laser detuning such that 
A(^) = ri(^~^ — 5^). In addition we can make some 
of the two-body terms equal by choosing V{£,) = 8fl£,. 
These two choices define a manifold (^-manifold) in pa- 
rameter space (A(^), V{(,)) along which the Hamiltonian 
Q becomes H{^) = H[n, A{^),V{i)]. Note that along 
the ^-manifold the Rydberg Hamiltonian and iJns are 
similar but not equal, i.e. i?HS = -ff(C) + Ai/. The 
difference operator AH contains binary interactions be- 
tween Rydberg atoms at distance 2a as well as three and 
four-body interaction terms that are all completely ab- 
sent in the Rydberg Hamiltonian and, thus, cannot be 
accounted for by an appropriate parametrization. 

To see how close nevertheless the coherent state Q 
resembles the true many-body ground-state \G{£,)) of 



H on the ^-manifold, we compare the numerically de- 
termined ground-state energy E = {G{£^)\H{£^)\G{£,)) , 
as well as the average fraction of excited atoms / — 
^~^Efe,m(G(OI'^fe,m|G(^)) with the energy and excita- 
tion fraction of the Rydberg gas calculated with |^). In 
addition, we compute the overlap |(G(^)|^)p between the 
exact ground state and hard squares coherent state. Fig. 
[2] shows the data for a 6 x 6 lattice with periodic bound- 
aries. For small ^, i.e. low Rydberg density, both wave 
functions yield identical energies and lattice occupancies. 
This agreement is also reflected in the overlap integral, 
which shows that the many-body ground state of the Ry- 
dberg gas in the low density regime is given exactly by 
a PEPS of the form Q. In the opposite limit, ^ — 00, 
the energies differ signiflcantly, due to the difference AH 
of quantum hard squares and Rydberg gas Hamiltonian. 
However, we observe reasonable agreement of the calcu- 
lated densities as well as in the overlap. This is because 
for large ^ the ground state energy of both the quan- 
tum hard-squares model as well as of the Rydberg gas 
is minimized by the configuration with maximal occupa- 
tion of the lattice with hard squares. For intermediate ^, 
however, both calculations give different results for the 
fraction of excited atoms and the overlap integral is small. 

The remarkable point is that the exact ground state 
|G(^)) of -ff(^) is still well approximated for any ^ by a 
coherent state of the form Q. To see this, we introduce 
an independent, variational parameter i], write the PEPS 
in terms of this parameter and optimize \ri) by minimiz- 
ing the energy functional {r]\H{(,)\r]) . To compute the 
functional, we express it in terms of expectation values 
of classical operators and numerically calculate them via 
the partition function Z^i using a transfer matrix method 
j27] . Comparing the exact ground-state energy with the 
variationally obtained one (Fig. [2]) we see excellent agree- 
ment for all ^ as well as a very good agreement for the 
fraction of excited atoms. The simple optimization pro- 
cedure also drastically improves the overlap, which is well 
above 0.5 for most values of ^ and never drops below 0.25 
which is remarkable for a many-body state in a Hilbert 
space of dimension 67022 |5Sj. This result suggests that 
the 2D Rydberg gas permits the preparation of a PEPS 
of the form Q and that the hard squares coherent state 
Q is well suited for carrying out a variational study of 
the phase diagram of the Rydberg lattice gas. 

To determine the phase diagram of the Rydberg gas 
in the whole parameter range (Q, A,V) we minimize 
{H),j = {ri\H\ri) with H given by eq. ([I]). The posi- 
tion at which the functional becomes minimal is denoted 
by rjmin- The system is in the disordered (liquid) phase, 
if ?7min < Ccrit, whcrc ^crit is the transition point of the 
hard squares model. For rj^-m > ^crit it is in the (anti- 
ferromagnetically) ordered phase with broken sub-lattice 
symmetry [c.f. Fig{ljb)]. 

The results of our variational approach are summarized 
in the phase diagram shown in Fig. [3] as a function of the 
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FIG. 3. (Color online) Variational phase diagram of the 2D 
Rydberg lattice gas with nearest-neighbor blockade on a 16 x 
16 square lattice with periodic boundaries as function of the 
scaled detuning 5 = A/l^ and Rabi frequency uj = Q./V . The 
dashed red line depicts the ^-manifold along which the data 
of Fig. [2] have been determined and the dashed black line the 
values where »)niiii = Ccrit for u} < 0.1. All other features are 
explained in the text. 



scaled detuning S — A/V and laser driving oj — fl/V. In 
the region w > wi w 0.1, the system exhibits two distinct 
phases, a liquid and a solid phase. At large negative de- 
tuning, the Rydberg gas favors the solid phase near half 
filling (/ < 0.5). The laser driving acts against the ten- 
dency of the atoms to order, such that if \S\ is reduced 
the system undergoes a first order phase transition from 
the ordered phase to the liquid state (solid black line in 
Fig. [3]). Near the phase boundary the energy functional 
{H)rj has a double well structure with the minima located 
on either side of ^crit- The transition is discontinuous, as 
the minima change their relative heights when passing 
through the transition (c.f. inset of Fig. [s] which ex- 
emplifies this behavior along the line A — B — C). For 
wi > w > aj2 ~ 0.07, we observe a qualitative change of 
the system's behavior. Close to the phase boundary, the 
energy functional still has two distinct minima, however 
both are located below ^crit , indicating a first order tran- 
sition from a low- to an high-density phase when reducing 
the detuning below \S\ « 0.4. At a; = W2 the transition 
line appears to terminate at a critical point. Here the two 
wells of the variational energy {H)jj merge into a single 
minimum. Below uj = uj2 our variational state jry) pre- 
dicts a continuous crossover from the liquid to the half 
filled state with increasing negative detuning. In the limit 
w — > we expect two classically ordered phases, one at 
quarter filling for > (5 > —0.5 and the staggered phase 
at half filling for —0.5 > S. The variational state \ri) can- 
not resolve these expected discontinuous transitions for 
u! ^ UJ2- We emphasize that the discussed features of the 
phase diagram are a consequence of the long-range (be- 
yond nearest neighbor) Rydberg interaction. For strict 



nearest neighbor interaction the variational calculation 
predicts a significantly shifted transition line and a first 
order phase transition from the empty lattice to a fully 
crystalline state starting at 5 = lj = (c.f. yellow line 
in Fig. |3|. We have checked that including interac- 
tion terms that go beyond next-nearest neighbors only 
marginally change the phase diagram. 

Our method for exploring the strongly correlated 
ground state of a 2D Rydberg gas appears to be fruit- 
ful also for other lattice geometries. Aside from studying 
the phase diagram it will be interesting to further ex- 
plore the resulting PEPS, e.g., with focus on their power 
to perform computational tasks [221 126| . 
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